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PREDICTION OF THE EFFECTS OF RROIELLER OPERATION CEN THE 
STATIC IDNGITEDINAL STABIL IT Y OF SINGLE-ENGINE 
TRACTOR MONOPLANES WITH FLAPS RETRACTED 
By Joseph Nell and. William C. Sleeman, Jr. 

I 

SUMMARY 


The effects of propeller operation on the static longitudinal 
stability of single-engine tractor monoplanes are analyzed, and a Bimple 
method is presented for computing power— on pitching-moment curves for 
flap— retracted flight conditions. The methods evolved are based on the 
results of powered-model wind-tunnel investigations of 28 model configu- 
rations. Correlation curves are presented from which the effects of 
power on the downwash over the tail and the stabilizer effectiveness can 
be rapidly predicted. The procedures developed enable prediction of 
power-on longitudinal stability characteristics that are generally in 
very good agreement with experiment. 


INTRODUCTION 


The prediction of the effects of propeller operation on the static 
longitudinal stability and control characteristics of single— engine 
tractor airplanes has been the object of many investigations. Successful 
methods have been developed for estimating the direct propeller forces 
and the effects of slipstream on the wing— fuselage characteristics 
(references 1 to 4) . Attempts to predict the complex changes In flow 
at the tall plane, however, have been somewhat less successful, primarily 
because many of the early researchers were hindered by insufficient 
experimental data for developing methods of proven general applicability. 

During the war yearB an appreciable amount of experimental data 
pertaining to power effects on static longitudinal stability was obtained. 
An analysis of these data suggested the possibilities of a semiempirical 
approach to the problem of determining the effects of power on the tail 
contribution to stability. This approach has been followed in the 
present paper nnfl a simple, rapid method for determining the effects of 
power on the tail contribution is presented. Use of the procedures 
developed permit the accurate prediction of power— on longitudinal 
stability and trim characteristics. No analysis has been male for the 
flap-deflected condition. 
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SYMBOIS 

lift coefficient 
pitching-moment coefficient 

average section pitching-moment coefficient ahont 
aerodynamic center for wing section immersed, in 
slipstream 

thrust coefficient 

Vpy^ 

thrust coefficient corresponiiing to power-off lift 
coefficient 

increment of thrust coefficient from power-off condition 
to a specified power condition 

airspeed, feet per second 

air density, slugs per cubic foot 

propeller disk area, square feet 

area of wing or tail, Bquare feet 

span of wing or tail, feet 

propeller-blade section chord, feet 

propeller diameter, feet 

propeller radius, feet 

radius to any propeller hlade element, feet 

wing mean aerodynamic chord, feet 

wing root chord at plane of symmetry, feet 

wing chord at break for wings having conqpos ite plan 
forms, feet 

wing chord at theoretical tip, feet 

wing chord at spanwise station O. 5 OK or 0.75R from airplane 
center line, feet 
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A 

wing aspect ratio 

X 

wing taper ratio (ct/ c r for vings having linear taper) 

z 

distance from, reference center of gravity to thrust line 
measured perpendicular to thrust line (positive 
when c.g. is above thrust line), feet 

h 

distance from reference center of gravity to propeller 
center line measured parallel to thrust line, feet 

h 

distance from reference center of gravity to elevator 
hinge line, feet 

h 

distance from elevator hinge line to thrust line 
measured perpendicular to thrust line (positive 
when elevator hinge line is above thrust line), feet 

a 

nngl« nf attack, radians unless otherwise denoted 

0 

propeller blade angle, degrees 


stabilizer setting with respect to thrust line (positive 
when trailing edge is down), degrees 

&e 

elevator setting with respect to chard line of Btabilizer 
(positive when trailing edge is down), degrees 

e 

effective angle of dcwnwash at horizontal tail, degrees 

e' 

increment of power-off downwash at horizontal tail from 
zero lift downwash, degrees 

^=0 

pcrwer—off downwash angle at zero lift 

Cyi 
^ * 

derivative of propeller normal— force coefficient with 
respect to angle of inclination of thrust- line in 
radians 

°r> 

v o 

value of Cyi for T c = 0 

Tjr 


S.F.F. 


abbreviation for propeller side— force factor 

(prf .2 ir sto ( p - P °-'?5E + 25 °) a §) 


f 

fr 


ratio of Cyt^ to 0y» 


*r 


ratio of Cyt^ for power-off value of T c to Cyt 
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Subscripts: 

T 

e 

t 

P 

P 

0 

V 

wf 

1 


empirical taper-ratio factor 

ratio of power-on stabilizer effectiveness 


to power— off stabilizer effectiveness 


&, 
Alc„ 


ratio of power—on elevator effectiveness ( — — 

V d5 e, 

to power-off elevator effectiveness 


change in a quantity due to power 


thrust line 
elevator 
horizontal tail 
propeller 
power on 
power off 
wing 

wing-fuselage combination 
immersed in slipstream 

BASIS OF ABAETSIS 


The method of computing power—on pitching moments which is 
outlined herein is based on the assumption that power-off (propeller- 
off or w indmilling ) pitching-moment and lift data are available for at 
least two stabilizer settings and with the tail off. The accuracy 
with which the effects of power on the tail contribution to stability 
can be predicted is dependent on the basic power-off data, and when 
possible these data should be obta i ned from wind-tunnel tests. 
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When pcwer—off wind— tun riel data are not available for use in preliminary 
design, the powe2>-off characteristics may he estimated hy use of references 5 
to 11. The wing mean aerodynamic chord and aerodynamic center may he found 
hy the method presented in reference 5- The lift— curve slope, angle of zero 
lift, and pitching-moment characteristics of the wing may he computed hy use 
of references 6 and 7* The effect of the fuselage on the wing— fuselage 
pitch in g moments may he determined hy Multhopp*s method “(references 8 
and 9) • A satisfactory approximation of the horizontal— tail effectiveness 
can he obtained when the isolated horizontal— tail effectiveness found hy the 
method of reference 10 is multiplied hy a factor of 0.9. 

The variation of power-off downwash with angle of attack computed hy 
use of the charts of reference 11 generally was found to agree fairly well 
with the variation of effective downwash with angle of attack obtained 
from wind-tunnel data when the computed downwash was multiplied hy a 
factor of 0.9 for all conditions for averaging downwash across the tail 
span instead of the factors obtained from figure 21 of reference 11. The 
absolute angle of downwash computed by use of the charts of reference 11 
had to he adjusted, of course, so that this angle would agree with the 
test downwash angle at zero lift. This adjustment was necessary since 
an appreciable amount of effective downwash was found to exist at zero 
lift chiefly as a direct result of local flow angularity at the tail 
caused hy the flow pattern over the rear of the fuselage. This downwash is 
often difficult to predict accurately. Neglecting the downwash at zero 
lift, however, will not affect the basic longitudinal stability or the 
estimated power effect® hut will alter only the trim characteristics. 

The experimental data upon which the results of this paper are based 
were obtained from wind-tunnel investigations of powered models of specific 
military airplanes. In figure 1 two viewB of each model and in table I 
the geometric characteristics of the configurations used are presented. 

Most of the models were tested in the Langley 7 ~ by 10— foot tunnel at an 
effective Keynolds number of approximately 1.6 x 10°. The power-off data 
were obtained with the propeller windmilling at a value of T c ~ -0.015 . 
Models 25 to 27 were tested in the Ames 7~ by 10— foot tunnel at an 
effective Beynolda number of approximately 2.0 x 10^. The basic power- 
off data used for these latter models were obtained with the propeller 
removed. 


METHOD OF ANAETSIS 


In the following discussion the individual component effects 
contributing to the over-all power-on static longitudinal stability are 
treated separately and approximate formulas are developed for estimating 
these effects. 
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Effect of Power on the Wing-Fuselage Characteristics 

Direct propeller effects .— The increment of lift coefficient 
contributed by the direct propeller thrust due to the inclination of the 
thrust line (the lift component of the normal force which is us ually small 
being neglected) is given by the following equation: 


AC Lp = T c sin arp (1) 


The increment of pitching-moment coefficient contributed by the 
propeller as a direct result of the thrust and the normal forces is given 
by the following equation, which was developed from equation (5) of 
reference 1: 




h - at_ 

it Cw c 


(f - fo)°y»^ 

r 


V 


( 2 ) 




where is the absolute angle of attach (radians) of wing from zero 

lift. Figures needed for use in equation (2) have been reproduced from 
reference 1 and are presented as figures 2 to 5. The term f — f r 
obtained from figure 2, is the difference in Gyt^Oyt^ 

and power-off conditions. It should be noted, however, that 
when power-off data are obtained with the propeller removed. 


L o* 

for power-on 


fp = ° 


term Oyt 




The 

obtained from figures 3 and 4, is the normal— force 


derivative; figures 3(a) and 4(a) are for low-speed propellers having 
thick, cambered blades; figures 3(h) and 4(b) are for high-speed 
propellers having thin, wide blades; plan-form curves of propellers on 
which figures 3 and 4 are based may be found in figure 4 of reference 1. 
The term de/da, the upwash factor, is obtained from figure 5* 


Slipstream effect on wing— fuselage characteristics .— The method 
most widely used for computing the increase in wing lift due to the 
slipstream is given by Smelt and Davies in reference 3. This method 
requires several successive approximations, however, to obtain final 
power-on lift coefficients when T c varies with Cj,. An approximate 
formula has been developed which is shorter than that of reference 3 
and which requires only a Bingle estimation to obtain the final 
value of AG^; thus an appreciable amount of computing time is saved. 

This equation is given by 



HACA TN No . 1722 


7 


A %r = 0.57T ca C l 


fw d£ 

O Off Sy- 


(3) 


where is the wing chord at epanwise station O.75R from airplane center 

line for wings "behind single rotating propellers or 0 . 50 R for wings "behind 
dual rotating propellers. 

Very close agreement is shown "between the values given "by equation (3) 
and two approximations computed "by the method of reference 3. (See fig. 6.) 
Somewhat less agreement is shown "between the values given "by equation (3) 
and test data (fig. 7) . The scatter Bhown can "be attributed to the 
idealized assumptions in the theory of reference 3 and the experimental 
inaccuracy of the test data. The effect of propeller tilt on AC]^. is 

small as shown "by reference 12, the data of model 2k, and other unpublished 
data and may be neglected for tilt angles up to at leaBt 5°. 

The effect of the slipstream on wing pitching moments is Hm»n in 
some cases, but it may be relatively large in others. This pitching-moment 
increment is obtained from equation (5) of reference 2 and is given as 
follows: 


°w 


AC™ = Cm -Ai + 



w 


where 

^ac 



average section pitching-moment coefficient about 
aerodynamic center for part of wing immersed in 
slipstream 

span of wing immersed in slipstream (taken as O.9D) 
average chord of wing immersed in slipstream 


rate of change of wing-fuselage pitching-moment 

coefficient with lift coefficient (propeller off) 
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Effect of Power on -the Ta ll Contribution to Stability 

Change in downwash angle due to power .— The downwash at the tail 
plane with the propeller removed is known to be chiefly dependent -upon 
the wing lift coefficient and the location of the tail with respect to 
the wing vortex system.. "When a propeller is added in front of the wing, 
many complex c hang es in flow occur which affect the downwash over the 
tail, but the chief effects are probably caused by the increase in 
wing lif t coefficient and altered wing span load distribution brought 
about by the passage of the propeller slipstream over the wing. Although 
.appreciable downwash may exist behind' an isolated propeller* at an 
angle of attack, large changes in the inclination of the thrust line 
(at constant wing angle of attack) -were found (reference 13 ) to cause 
practically no change in effective downwash at the tail when a wing 
was located between the propeller and the tail. "With the foregoing 
discussion as a basis the following simplified semiempirical approach 
was used to derive a parameter with which to correlate experimental 
downwash changes due to power. 

Downwash angles were emigrated for several models for which 
extensive constant thrust data were available. When Ae was plotted 
against T c at various angles of attack, Ae was found to be a 
function of both T c and a. This relationship seemed logical for 
the increase in wing lift with power, and any resultant span— load 
changes would also depend on T c and <x> however, e* was believed 
to be a better factor than a for use in the correlation inasmuch 
as e’ usually varies linearly with a up to fairly high lift 
coefficients and also depends on tail location. The aBBungition was 
made that a tail well out of the power— off maximum downwash field 
would also be favorably located in the power-on downwaBh field for 
configurations within the range of geometry of the Aaodels presented. 

Model 17, which has an untapered wing, showed a considerably larger 
increase in e with the application of power than either model 13 or 
model l6 , which had identical tail and fuselage geometry but rather 
highly tapered wing plan forms. The power-off downwash angles were 
considerably less for the model with the untapered plan form, but the 
downwash for all three wing plan forms could be accurately computed 
from the charts of reference 11. With power, however, the downwash 
angles for models 13 , l6, and 17 were much closer to the same value. 

According to lifting— line theory the downwash behind a wing of 
arbitrary plan farm in a uniform air stream (at a given tail location 
and lif t coefficient) depends only on the span load distribution along 
the wing. Wing taper ratio has a significant bearing on the span load 
distribution and hence the downwash at the tail, the downwash increasing 
with wing taper. The data of reference 14 show that the slipstream 
alters the span load distribution by increasing the load in g over the 
inboard part of the wing. Since the greatest change in downwash with 
taper occurs for values of 1 near unity, a wing of rectangular plan form 
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might he ejected to be more susceptible to this induced taper effect 
and show the largest change of downwash angle with power. The experi- 
mental results of models 13 and 17 agree with the foregoing assumptions; 
however, further substantiation would be desirable. 

An empirical factor F (fig. 8) which is a function of wing taper 
ratio was derived from the data of models 13, 16, 17, and other models 
with s imil ar tail geometry to account for induced taper effects. The 
taper ratio for wings of composite plan form may be satisfactorily 
estimated by use of an equivalent “root chord c r * , aB ■ shown in the sketch 
in figure 8. 

A plot of the parameter (AT c )e’F against the experimental Ae 

obtained from stabilizer and tail-off wind-tunnel data for 28 model 
configurations is shown in figure 9 (a) . The dash— line curves in the 
figure indicate the approximate accuracy of determining downwash # angles 
from complete-model wind-tunnel data. The correlation of test points 
indicates that the parameters selected account for the first-order 
effects of power rather well. The solid— line curve indicates the 
suggested curve for use in design and, is reproduced in figure 9(b) with- 
out experimental test data. 

f 

Rote that figure 9 (t) indicates no change in Ae attributable to the 
tilt of the propeller thrust axis. The data of model 2 h and other unpub- 
lished data show that changes .in Ae with tilt are small and rather 
inconsistent and the effect of tilt (at least for tilt angles up to 50) 
can be satisfactorily estimated from, considerations of direct thrust 
effects and changes in stabilizer effectiveness. 

Change in stabilizer effectiveness with power .— The slipstream is 
considerably distorted in the region of the horizontal tail, and 
idealized theoretical methods which assume a cylindrical slipstream at 
the tail do not always produce a satisfactory estimate of the change 
in dynamic pressure. at the tail associated with the application of power. 

As was true for the downwash correlation, a semi empirical approach was 
followed to derive a method for predicting the change in stabilizer 
effectiveness with power. 

The ratio of power-on to power— off stabilizer effectiveness 
was assumed to be directly proportional to AT C and, the ratio of the 
propeller diameter to tail span. A maximum value of E^. was also assumed 
to be attained for the tail located on the thrust line with a linear 
decline in E-^ occurring until a value of unity is reached for a tail 
location 1 propeller diameter above or below the thrust line. Experimental 
points were plotted against the paramecers suggested by the foregoing 
assumptions (fig. 10) and the following relationship was obtained; 



B -|- = 1.0 + 2.1 


( 5 ) 
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The dash-line curves in figure 10 indicate the appr oxima te limits of 
accuracy of determining B t from wind— tunnel test data, and essenti ally 

all the test points fall within these limits. 

Change in elevator effectiveness with power .— Inasmuch as most 
airplanes utilize the elevator for longitudinal control, it is apparent 

that the increase in elevator effectiveness ~r^ with power will influence 

d 6 e 

the determination of the final power-on stability and trim c har acter! sties 
of the airplane. An analysis was made to determine the possibility of 
correlating B e in a manner similar to the foregoing correlation of B t . 
The results of this analysis were less consistent than the results obt ain ed 
with the correlation of B^ , probably for the most part because of 
appreciable scale effects on some of the model elevator data the 
reduced accuracy possible in setting the elevator-deflection angles. For 
full-span elevators when estimated power—on elevator data are desired, 
however, it may be assumed in most instances that E 0 = B^.. 


Computation of Power-on Lift and Pitching-Lament Coefficients 

Power-on wing-fuselage lift coefficient .— The final power-on wing- 
fuselage lift coefficient is given by 


(U = (M, + ( a %) 2 + 


(6) 


In order to arrive at a value of 


from equation '( 6 ), the 


following procedure is recommended for conditions where T„ varies 

v 

with C^: The increment of the wing lift coefficient due to power is 

first evaluated by equation ( 3 ). The second-approximation of the Increment 
of lift coefficient contributed by the propeller ^ACLp ^ 2 iB obtained by 

computing a first approximation ^ACpp^ by equation (l) with values 
of T c correspo ndi n g to + AC]^j the second approximation is 

then obtained by use of equation (l) with values of T c corresponding 
to (OLvz) + + ^ACppj^. The second approximation will usually 

give a value of such that further approximation will be 


unnecessary . 
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P over-on wing— fuselage pitching-moment ooef f lclents Tie final 
power— on wing-fuselage pitching-moment coefficient is given by 


( < W) TJ " ( Cm wf) 0 + ACm P + ACm v 


( 7 ) 


The terms ACmp and AC^ are found from equations (2) anfl (4), with 
values of T c based on as given hy equation (6) . 


Pover-on tail pitching-moment coefficient.— The computation of the 
power-on tail pitching-moment coefficient merely consists of adding the 
increments produced hy the altered downvash at the tail and increased 
tail effectiveness to (^m^j * this coefficient is given hy 


(Si = Et (S ) 0 As [ Bt © 


( 8 ) 


Power-on complete-model pitching-moment and lift coefficients .— 
The final power-on complete-model pitching-moment coefficient iB given 
hy adding equation (7) and equation (8) as follows: 


Inasmuch as 


% = (^*f)p + (S) p 



( 9 ) 


( 10 ) 


the final power-on complete-model lift coefficient is given hy 
adding equation (6) and equation (10) as follows: 


% ' (V) p ' + (S) P 


(ii) 
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Illustrative Example 

A detailed step-b>y— step procedure for computing power— on lift 
coefficients and pit ch i n g-moment coefficients for model 21 is presented 
in table U. The sample calculations in table XX illustrate the use 
of the equations presented in this paper and give the pertinent 
c ons tants and column headings in a convenient form for general application 
to design. The data from which these estimations were made and the final 
computed power-on characteristics are presented in figures 11 and 12 . 

Model 21 was chosen as an example because* although the individual 
component effects of power were not small* the design variables were 
such that the adverse effects were counteracted by the favorable effects 
ami thus a very ama .1 1 over— ell change in stability with power resulted. 
Calculations for models 13 and 15 were also made to show that the change 
of power effects attributable to raising the horizontal tail and 
increasing Its area can be accurately predicted. The basic data and 
estimations of power effects for these models are presented in figures 11 * 
13 * 'and Ik. 

The computed results for al 1 three models show very good agreement 
with the test results. 


DISCUSSION 


The range of the most pertinent geometric variables for the models 
used in this paper are presented in table ill. The correlation curves 
of figures 6 * 9 * and 10 are believed to be valid at least between the 
limits given in table TTT . No data on powered models with appreciable 
wing sweep, were available* consequently* the effect of sweep could not 
be included in the correlating parameters. 

Wind-tunnel data on personal-type airplanes were not available for 
use in the correlations* and the applicability of the curves in 
figures 9 and 10 to this type of design is dependent on a number of 
factors. Although the models used in the present correlation represent 
highly powered fighter-type airplanes* the correlation curves were 
found to be valid for medium power conditions on the fighter-type 
airplane models also. An estimate of the variation of T c with C^ 

for several typical single— engine personal-type airplanes showed that the 
thr ust coefficients for maximum rated power for these airplanes fell in 
the range of thr ust coefficients for the medium power conditions on the 
fighter-type airplanes. 

The range of wing vertical positions relative to the slipstream, 

- and the ratio of the slipstream diameter to the wing span might be 
expected to be considerably different for -military and personal— type 
airplanes* and these differences could have a significant bearing on 
the magni tude of the power effects. 
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All models presented herein have a wing location that placed the 
■wing well within the slipstream. The increment of wing lift due to the 
slipstream derived from the data of reference 15 does not vary 
appreciably with wing height for wing positions 0.3 propeller diameter 
above and below the thrust line when the propeller is more than 0.3 root 
chord ahead of the wing leading edge. The range of wing vertical positions 
for the models presented herein is O.165 and O.I76 propeller diameters 
above and below the thrust line, respectively, and the propellers are 
more than 0.3 root chord ahead of the wing leading edge; thus the range 
falls within the l im its of wing and propeller locations where computed 
values of ACt would be expected to -be valid. 


Generally, the diameter of the propeller relative to the wing span is 
smaller for personal— type airplanes than for fighter-type airplanes. 

Model 2 k has a relative propeller diameter approximately the same as for 
the personal-type airplanes considered, but the other models used in the 
correlations had larger relative propeller diameters. Bo definite 
conclusions can be made regarding the effect of relative propeller size 
because of insufficient data. In most instances, the methods outlined 
in the present paper should be satisfactory for computing the first- 
order effects of propeller operation on personal-type single-engine 
tractor airplanes. 


OPTIMUM 33 ESIGB CONSIDERATIONS 


The design configuration usually considered optimum when satisfactory 
handling qualities of airplanes are considered is that which exhibits no 
change in longitudinal stability characteristics upon the application of 
power. Many design parameters affect the longitudinal stability both 
adversely and favorably, and defining a definite method by which to design 
an airplane with no power effects iB difficult. Often considerations 
other than aerodynamic determine the final geometry of a design. In view 
of this fact and the rapidity with which the power effects of a specific 
configuration can be computed by the method of the present paper, each 
proposed design should be examined for power effects, and an optimum 
configuration {minimum power effects) should be attained by a process of 
rational modification to the original design. 


Langley Aeronautical Laboratory 

National. Advisory Committee for Aeronautics 
Langley Field, 7 a., July 13 , 19^5 
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_ TABLE HI.- RANGE OF GEOMETRIC CHARACTERISTICS 
OF MODELS INCLUDED IN CORRELATION 


Geometric parameter 


Range 



Wing aspect ratio 


Wing taper ratio 


Propeller diameter 


Wing span 


Height' of tail above thrust line 


Propeller diameter 


Tail span 


Wing span 


Tail length 


Mean aerodynamic chord 


Height of thrust line above vlng root chord 


opeller diameter 


Distance of propeller ahead of vlng root chord 


Root chord 
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NACA TN No . 1722 
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Figure 1.- Continued. 
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Figure 1.- Concluded. 
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0 ID £0 30 bo 50 to 

p at 0.75 R, deg 

(a) Hamilton Standard 3155-6 propeller with 0.164 -diameter 
s pinn er; S JJ. m 80.7. 



0 ID 2 D 30 bo 50 Go 

P at 0.76 R, deg 

(b) NACA propeller 10-3062-045 with 0.164 -diameter 
Epinnar; S.F.F. ■ 131.6. 


Figure 3.- 


Variation of C v i 
x tJt 
¥ 0 


with blade angle. 






6 blade, dual rotation 
,2 b l ade, single rotation 

3 blade, single rotation 

4 blade, single rotation 
6 blade, single rotation 
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HA.CA r TN No. 1722 
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Figure 7.- Comparison of AC T computed by equation 3 and obtained from 


test data. AC 
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Effective thrust coefficient 


1.0 


.8 

o 
Eh 

.6 
.4 
•2 
0 

Figure 11.- 
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Lift coefficient, 

Variation of constant power thrust coefficients with lift coefficient for models 13, 

15, and 21. 
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Lift coefficient, 

(b) Comparison of computed and experimental data for constant power condition. 


Figure 12.- Continued, 




L.2 



HACA TN No . 1722 


a 

o 


.1 


» 


-8 


Fii 


Angle of attack, » , deg Pitching-moment coefficient, 










Shift in neutral point location due to power, An D , percent 
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(c) Comparison of computed and experimental effects of power on the neutral point 

and tail-off aerodynamic center location. 


Figure 13.- Concluded. 
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